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椭圆 型 方程 的 一 种 新 型 混合 有 限 元 格式 * 


We, FEFE, AR 
(西安 交通 大 学 理学 院 ， 西 安 710049) 


摘 要: 基于 实际 问题 对 通 量 较 低 的 正则 性 要 求 ， 本 文 建立 了 椭圆 型 方程 的 一 种 新 型 混合 变 分 形式 . 在 鞍 
点 问题 框架 下 ， 通 过 验证 LBB 条 件 ， 证 明了 该 混合 形式 解 的 存在 唯一 性 .由 于 压力 空间 不 再 是 
传统 的 H (div)， 而 是 平方 可 积 空间 ， 因 此 混合 元 的 选取 变 得 简单 容易 .同时 本 文 给 出 了 相应 的 
有 限 元 逼近 形式 ， 并 对 于 由 分 片 常数 速度 元 和 分 片 线性 压力 元 构成 的 协调 有 限 元 对 遍 -已 ， 通 
过 验证 速度 投影 算 子 的 有 界 性 ， 证 明了 有 限 元 解 的 存在 唯一 性 ， 以 及 有 限 元 允 近 在 某 种 意义 下 是 
最 优 的 . 最 后 给 出 了 数值 算 例 ， 验 证 了 方法 的 有 效 性 和 理论 分 析 的 正确 性 .此 方法 还 可 以 通过 增 
加 简单 的 稳定 项 使 用 最 为 常用 的 最 低 阶 等 阶 有 限 元 . 

关键 词 : 椭圆 型 方程 ， 混 合 变 分 形式 ， 协 调 有 限 元 对 ; LBB% 

分 类 号 : AMS(2000) 65N12; 65N30 ”中 图 分 类 号 : 0175.25; O241.82; O242.21 ”文献 标识 码 : A 


1 引言 
NÆR? (d = 2,3) 空间 中 具有 Lipschitz 边界 (T = 00) 的 有 界 区 域 . 我 们 考虑 如 下 的 椭圆 型 
方程 
—Ap=f, 在 9 内， (1) 
p=0, ÆT E. (2) 


该 方程 刻画 了 很 多 的 物理 和 力学 现象 ， 例 如 : 具有 固定 边界 的 弹性 膜 在 强加 强度 为 了 的 外 力 时 
的 平衡 状态 ， 当 f = 0 时 ， 一 些 物 理 量 (化 学 浓度 、 温 度 分 布 等 ) 的 平衡 问题 等 由 ,对 该 方程 的 
有 限 元 逼近 方法 和 为 了 降低 解 的 正则 性 要 求 而 引入 的 混合 变 分 形式 及 其 有 限 元 逼近 已 经 有 许多 
ALA, 

假定 fe L2(0)， 引 入 通 量 u = 一 Vp， 方程 (1)-(2) 传统 的 混合 变 分 形式 为 : Kuev, pe 
W， 使 得 


(u,v) -—(V-v,p)=0, VveEV, (3) 

(V-u,w)=(f,w), Vwew, (4) 
其 中 

V = H(div, 9), W = L?(Q). (5) 
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令 Kh 是 的 一 个 剖 分 ， 选 取 适 当 的 Vi C VAW, C W， 则 (3)-(4) 相应 的 有 限 元 逼近 
为 : Ke up, E Vh, Ph € Wa> 使 得 


(up, v) — (V: v, ph) =0, Vv Ee Va, (6) 

(V -up w) = (f, w), Vw EW. (7) 
众所周知 ， 如 果 有 限 元 对 (Vh, Wah) 满足 离散 LBB 条 件 
(V e Vh, Wh) 


inf sup ———~— 

wnEWh v eVa llvallvl|wallw 

SE ERRERA, 则 (6),(7) 存在 唯一 解 ， 且 逼近 是 最 优 的 . 满足 要 求 的 有 限 元 对 有 很 多 ， 比 
如 RT 元 、BDM 元 等 等 2.3] . 

我 们 注意 到 通 量 u = 一 Vp 并 不 需要 (3)-(4) 中 如 此 高 的 正则 性 ， 事 实 上 仅 需 要 具有 LE 
性 ,基于 此 ， 我 们 在 下 一 节 中 我 们 给 出 (1)-(2) 的 一 个 新 的 混合 变 分 形式 ， 并 证 明 其 解 的 存在 唯 
一 性 . 在 第 3 节 中 我 们 给 出 相应 的 混合 有 限 元 逼近 ， 并 对 由 分 片 常数 速度 元 和 分 片 线性 压力 元 
构成 的 协调 有 限 元 对 证 明 解 的 存在 唯一 性 . 另外， 我 们 还 证 明 有 限 元 解 的 逼近 在 某 种 意义 下 是 
最 优 的 . 最 后 ， 在 第 4 节 中 给 出 一 个 数值 算 例 ， 验 证 第 3 节 中 给 出 的 有 限 元 逼近 的 有 效 性 以 及 
理论 分 析 的 正确 性 . 


>ß>0, (8) 


2 一 个 新 的 混合 变 分 形式 
由 分 部 积分 公式 和 p|r = 0， 可 知 (Y:v,p) = 一 (Vv, Vp). 因此 (3)-(4) 可 以 改写 为 


(WVv)+(v, Vp)=0, vveV, (9) 
—(u, Vu) = (f, w), Vwe Wi, (10) 
其 中 
Vi = (1(0)", Wi = (0), (11) 
并 赋予 通常 的 范 数 . 
定义 
a(u,v)=(u,v), Vu,ve Vi, (12) 
b(v,p) =—(v, Vp), VvEVi, pEW, (13) 
根据 (9)-(10)， 方 程 (1)-(2) 的 一 个 新 的 变 分 形式 为 : 求 ue V1, pew, te 
a(u,v)—b(v,p)=0, VveVi, (14) 
b(u,w) =(f,w), Vwe Wy. (15) 


显然 ， 这 是 一 个 鞍点 问题 . 针对 此 问题 ， 我 们 在 下 面 的 两 个 引 理 中 给 出 一 些 相 关 性 质 . 
引 理 1 双 线 性 形式 a(:) 和 1(.,) 是 连续 的 ， 且 a(-) 是 椭圆 的 ， 即 


a(u, v) < ujlv, llvlly,, Vu,v E Vi, (16) 
b(v,p) < lvllvillpllw,, VveVi, pew, (17) 


a(u,u) = uy Vue Vi. (18) 
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证 明 是 直接 的 ， 我 们 这 里 略 去 . 
引 理 2 ” 双 线 性 形式 以 .,-) 满足 LBB 条 件 ， 即 存在 常数 B! > 0， 使 得 
in pe 
wEW!1 N: llvllv, \|wl|w, 
证 明 对 任意 的 we WwW, ， 取 v = -Vw， 我 们 有 veEVi， 以 及 
~—(v, Vw) _ 1Ywlv， 
Iv lv, lwll, wll. 


根据 全 范 lwla: 和 半 范 [Vwi 的 等 价 性 ， 即 


> By. (19) 





all Ywl: < lwl < call Vwan we Wi, 
我 们 得 到 
—(v, Vw) x 1 
livllvi Iwllw, 一 ce 
ae + - 即 得 结 
结合 合 引 理 1 Ga 2， 我 们 从 经 典 的 鞍点 问题 抽象 理论 &,3 可 知 如 下 定理 成 立 . 
定理 1 变 分 问题 (14)-(15) 存在 唯一 解 (up) € Vi x Wi. 
注 1 (i) 车 fe HQ), BAe (1)-(2) 的 唯一 弱 解 满足 


pe 到 (0O)，u:= -vpe(L2(D)) 


这 些 正 则 性 和 定理 1 的 正则 性 是 一 致 的 ， 这 使 得 我 们 的 变 分 形式 比 (3)-(4) 更 加 自然 
(i) 进一步， 如果 fe Z2(Q)， 相 似 的 讨论 可 知 


(u, p) € (HI (OQ) x (HO) (Ho (2)) 


3 ARTELE 
在 选取 了 两 个 协调 有 限 元 空间 Vin C V1 和 Win C Wi 之 后 ， 我 们 给 出 (14)-(15) 式 的 有 限 
元 逼近 形式 为 : Run E Vin, ph © Win， 满 足 
alUn, Vh) — (Va, pr) = 0, VvVaE Via, (20) 
b(un, wn) =(f,wr), Vu E Win. (21) 


X T Wa EK IB UE i A AS FF EEE, 我 们 要 求 选取 的 有 限 元 空间 对 满足 形 如 (8) 的 离 
散 LBB 条 件 ， 即 存在 不 依赖 于 的 常数 6。， 使 得 
—(va, Vwn) 
f eae M 
wr€Win we [Valv eal, 
本 节 我 们 仅 考 虑 二 维 情 形 ， 令 Ki 为 0 的 三 角 剖 分 ， 并 在 Ciarlet 意义 下 是 正则 的 是， 
WR (Vin, Win) SARIGI P — Pi» 具体 定义 如 下 


Vin = {va = (vin, van) € Vi: vin € Po(K), Y K € Kn,i= 1,2}, (23) 


> Bo > 0. (22) 


Wi, = {w E€ C°(O) (Wi: w€ P (K), VK € Kn}. (24) 
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引 理 3 ”有限 元 对 P? -已 满足 离散 LBB 条 件 (22). 
证 明 我们 借鉴 Fortin 的 思想 25 来 验证 条 件 (22)， 仅 需要 构造 有 界 算 子 In : Vi 一 Vin» 
使 得 对 任意 的 ve Vo A 


b(v — Hav, wh) = 一 (V — IIhv, Vw) = 0, Y wh E Win- (25) 


& wh 一 CQKZ 十 OK 对 任意 的 K € Kn; 其 中 ax, Bk E RAMs K 上 的 任意 常数 . (25) 式 
等 价 于 
(v — Iv, (ax, bK)" ) g =9, 
即 
ax | ((v)' — (IInv)')dady + px f ((v)? — (IInv)?)dady = 0. 
K K 


Har, Bk € 五 的 任意 性 可 得 
Ilnv |k = xf vdrdy, VK EK, (26) 
| 天 | Jx 
从 而 由 Halder 不 等 式 得 到 


1/2 
aviv, = { 5 [avleclliceucys} 
KeK, 


-{ 5 a (( Sow) + ( f deas) 


KEK, 
1/2 
<{ I f tmPasay+ f loaPardy)}” = vliv: 
KEK, K K 
这 说 明 算 子 I 是 有 界 的 . 
Soh, WAER, RAJTA, Ty, 满足 插值 误差 
lv=-Invllv <chllvilun@ye, Vv € (A*(Q))%, (27) 


FE e KRF h All v. 
定理 2 对 于 有 限 元 对 P?2 — Pi, 混合 有 限 元 变 分 形式 (20)-(21) 存 在 唯一 解 (un,pn) E 
Vin x Win， 并 满足 


Ju—unllv, + lp paw < e( iaf, la- valiv + inf, lp- wnlw), Cg) 
Ht (u, p) € Vi x WW 为 (14)-(15) 的 唯一 解 。 另外 有 


ju — ually, + lip- paliw, < ch(llullencowya + lella), (29) 
进一步 地 
lp — pallz2cay < ch? (lula + lpllazc)- (30) 
证 明 “从 鞍点 问题 的 一 般 理 论 &,3， 可 得 到 (28),(29)， 我 们 只 需 征 明 (30). 为 此 定义 双 线 性 
形式 
Blu, p; v, w) = a(u, v) — b(v, p) + b(u, w). 
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那么 ，(14)-(15), (20)-(21) 可 以 分 别 改写 为 
Bu,p;v,w) =(f,w), Bun, Ph; Va, wh) = (f, wa). 
从 而 有 正 交 性 公式 
B(u— Un, P— PhiVn Wwa) =0, V (vas wn) € Vin X Win. (31) 


AT WEA (30) sh, BRA AS Layton-Toshiba 77 7K'671, se XIHA: 求 (6, 和 © Vi x 
Wi, 满足 
B(v,w3;¢,A) = (we), V(v,w) € Vi x Wy, (32) 


这 里 e = Pap- pr, Pr 为 通常 的 投影 算 子 ， 满 足 
lw- Prwllza + Allw — Palla < ch? jullie, Yw e RAW. (33) 
显然 ， 由 第 2 节 的 讨论 可 知 ，(32) 存在 唯一 解 ， 并 且 满 足 正则 性 
[Oller 十 人 as < elfellre. (34) 
令 I= Inu- un， 并 在 (32) 中 取 (v,w) =(n,e), RATE 
el = B(n, e; $, A), 
另 一 方面 ， 在 (31) PR (va, wa) = (nd, Pad)» WE 
B(u -Iru +n, p— Pap +e; Und—¢+¢, PrA-A+A) =O. 
结合 (29) 和 (32)-(34) 式 ， 并 在 (32) RPE (v, w) = (u — Hau, p — Pap). RAGS 
lele = B(u— Inu, p — Pap; Und — 6, Prd — A) 
+B(u — Iau, p — Php; $, à) + Bln, e; IInp — Q, PA — 2) 
< ch? (hull + llpllaz) (oll + IAllz2) 


+(p — Pap, e) + ch? (lull: +llplas)(llglla + |All) 


A 


IA 


ch? (julm: + llpllz2)llellz2, 


从 而 有 
llellz2 < cp? (lull + llpllz2)- 


进一步 由 三 角 不 等 式 ， 即 证 (30). 

注 2 从 证 明 过 程 可 以 看 出 ， 该 方法 可 以 直接 推广 到 三 维 情形 ， 而 且 对 于 二 维 的 四 边 形 剖 分 
和 三 维 的 六 面体 剖 分 也 有 类 似 的 结论 。 此外， 该 方法 可 以 自然 推广 到 一 般 的 椭圆 型 偏 微分 方 
程 ， 并 可 以 考虑 高 阶 的 稳定 化 有 限 元 格式 . 
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4 算 例 





在 本 节 中 ， 我 们 给 出 一 个 简单 的 算 例 来 验证 上 一 节 中 给 出 的 混合 有 限 元 格式 的 有 效 性 
以 及 理论 分 析 的 正确 性 . 在 原始 方程 (1)-(2) 中 ， 选 取 区 域 为 单位 正方 形 [0,1J?， 压 力 真 解 
Wp = sin(zrz) sin(rry)， 外 力 取 卫 = —27? sin(zrz) sin(rry)， 显 然 ， 压 力 满足 零 值 边界 条 件 . 

为 了 使 用 有 限 元 方法 求解 ， 对 区 域 进 行 一 致 网 格 前 分 : 分 别 沿 z 和 8Y 方 向 1 等 分 ， 然 后 对 
每 个 四 边 形 沿 同一 方向 对 角 线 连 接 ， 便 得 到 需要 的 三 角 剖 分 ， 对 于 不 同 的 网 格 ， 采 用 有 限 元 空 
间 (23)-(24) 求解 混合 有 限 元 格式 (20)-(21). 表 1 给 出 了 混合 有 限 元 逼近 解 在 不 同 范 数 下 的 误差 
和 收敛 阶 数 ， 


Rl: 不 同 尺度 网 格 下 的 速度 和 压力 逼近 误差 及 收 伊 阶 








M lp 一 pn 2 收敛 阶 llp-Prll g: 收敛 阶 lu 一 azllrz 收敛 阶 
lipi g2 Holl z1 lullg2 
10 2.787 E-2 / 1.586 E-1 / 1.624 E-1 / 


20 7.031 E-3 1.987 7.958 E-2 0.995 8.153 E-2 0.994 
40 1.762 E-3 1.997 3.983 E-2 0.999 4.081 E-2 0.998 
80 4.407 E-4 1.999 1.991 E-2 1.000 2.041 E-2 1.000 


160 1.101 E-4 2.001 9.960 E-3 0.999 1.020 E-2 1.001 





REL BTU, WEA REE AH FI A A ARRETA A RHEE J TR id 
ARTN ERER, MAENE L? A H 范 数 分 别 达 到 了 理论 上 的 2 阶 和 1 阶 收敛 精度 ， 同 时 
速度 也 达到 了 理论 上 的 了 ?2 范 数 意 义 下 的 1 阶 收敛 精度 . 
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A New Mixed Finite Element Scheme for Elliptic Equations 
SHI Feng, YU Jia-ping, LI Kai-tai 


(School of Science, Xi’an Jiaotong University, Xi’an 710049) 


Abstract: In this paper, a new mixed variational formulation to the elliptic equation is given based 
on the less regularity of flux in practical problems, and the existence and uniqueness of solution to this 
formulation is shown under the framework of the saddle point problem by verifying the LBB condition. 
Since in the new mixed variational formulation, the pressure belongs to the square integrable space 
instead of the classical H (div), the choices of finite element pairs become simple and easy. Based on 
this new formulation, the corresponding conforming finite element approximation is addressed, and the 
existence and uniqueness of finite element solutions are obtained for PÊ — P, pair consisting of piecewise 
constant element for velocity and piecewise linear element for pressure, by proving the boundedness of 
projection operator for velocity. Moreover, it’s proven that the finite element approximation is optimal 
in some sense. Finally, some numerical experiments verify the effectiveness and theoretical results of 
our method. This method can be extended to utilize the lowest equal order finite element pair by 
adding some simple stabilization terms. 

Keywords: elliptic equation; mixed variational formulation; conforming finite element pair; LBB 
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